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An overview is given of analyses in DIS at HERA which confront the predictions of power corrections with measured 
data. These include mean values and distributions of 2-jet as well as 3-jet event shape variables and jet rates. 



1. Introduction 

\q ■ 

Event shape variables are known to exhibit rather large hadronisation effects, not treated by purely perturbative 
calculations, ft has been shown Q] that the size of these corrections varies as (A/Q), with Q being the hard scale 
of the process, for most of the event shape observables. This article focuses on power corrections in the approach 
initiated by Dokshitzer and Webber Q], since these are the most complete theory and often used in comparisons 
K-5 , to data. There is a wealth of studies performed in e + e~ annihilation which gives support to this concept, cf. to 
' contribution R002 in these proceedings. It is interesting to extend these studies to deep-inelastic scattering (DIS), 
in order to investigate the universality of the ansatz and to check for possible modifications of the hadronisation 
I i process due to thepresence of a proton remnant. After power corrections to mean values of event shapes variables 
| became available |3j, the HI and later also the ZEUS Collaboration published analyses which corroborated this 
£^ ' concept 0,0. A comprehensive review on the subject of event shapes can be found in 0|. 

When comparing the ep scattering to e + e~ annihilation, one moves from an s-channel to a t-channel exchange, 
with the negative four momentum transferred Q corresponding to the center of mass energy y/s. At HERA a large 
range of the scale is available in a single experiment, typically stretching 5 < Q < 115 GeV. The theoretical treatment 
is complicated by initial state singularities related to the incoming proton, which are absorbed in the parton density 
functions (pdf). In order to reject the proton dissociation part of an event (the remnant), the event shapes are 
D , calculated in the Breit frame of reference, where the separation between particles from the hard scattering and the 
O ,' remnant is clearest. The boost to the Breit frame is determined by the event kinematics and the boosted particles 
are separated into hemispheres with pseudorapidity r\ < (current hemisphere) and n > (remnant hemisphere). 
Usually only particles of the current hemisphere enter the event shape definition, where this hemisphere resembles 
. \ to some extend one half of an e + e~ event. In order to ensure infrared safety at all orders a minimal energy in the 
' current hemisphere is applied as part of the observable's definition, e.g. Ech > Q/10. 

2. Event Shape Variables 



The observables covered in the following may be separated into three categories: 2-jet event shapes, 3-jet event 
shapes and jet rates. The term "n-jet" denotes that at least n particles in the final state are needed for a non 
vanishing value of the variable. Note that particles from the proton remnant is not included in n. Common examples 
for 2-jet event shapes are thrust, jet broadening, jet mass and the C-parameter. For this class of observables the 
most advanced theory predictions are available, therefore many analyses concentrate on those. 3-jet variables which 
have been proposed are the out-of-plane momentum and the azimuthal correlation. While the theoretical predictions 
are yet less complete for these variables, additional insights to power corrections are possible due to sensitivity to 
hadronisation from a gluon and applications to hadron-hadron collisions. 

Closely related to event shapes are jet rates which make use of a jet clustering algorithm, such as the kt or the 
JADE algorithm. Here the algorithm is applied on particles on both hemispheres and the n-jet rate is defined as 
the maximum value of the cut off parameter y cnt where the event is clustered to n + 1 jets (+1 denotes the proton 
remnant). 
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A particular distinction made for 2-jet event shapes is only possible for DIS: between ones making use of the 
exchanged boson axis (mostly from a virtual photon) in their definition and those without that reference axis, as 
for the definitions used in e + e~ annihilation. Namely for thrust and jet broadening both variants are investigated, 
whereby the explicit use of the boson direction implies sensitivity to radiation into the remnant hemisphere through 
recoil effects on the current quark Q. 

The thrust variable r with respect to the boson axis is defined as 



and the thrust variable tc 



t = 1 - T with T = ^ Pz J h } , (1) 

T,h \Ph\ 



tc = 1 - T c with T c = max — " — , (2) 

n -r T,h\Ph\ 



The squared Jet Mass is calculated as 



where the direction n*T maximises the sum of the longitudinal momenta of all particles in the current hemisphere 
along this axis. 

The Jet Broadening is defined as 

r> Y^h \Pt,h\ / Q \ 

9 (2EJAI) 2 ■ u 

In the following the symbol po is used, which indicates that in the above definition the hadrons are treated as massless, 
replacing the energy Eh by the modulus of the 3-momcntum \ph\- This adjustment is made since the theoretical 
predictions assume the partons to be massless. Mass effects can be huge, especially for the jet mass observable. 
The C-Parameter 



3Eaj,' \ph\\Ph'\sm 



2 , 



c 2 ' ThApW ' (5) 

where Ohh' is the angle between particles h and h' . Note that in all the definitions the momenta are defined in the 
Breit frame and the sums extend over all particles in the current hemisphere. 
The out-of-event plane momentum is defined as 

^out=EV° ut i- 

h 

Results are presented in terms of K ou ^/Q. Here p° ut is the out-of-plane momentum of the hadron h with the event 
plane defined to be formed by the proton momentum P in the Breit frame and the unit vector n which enters the 
definition of thrust major: 

T M = max — y \ph-n\, n-P = (7) 

n 

To avoid measurements in the beam region, the sum indicated by ^2 h ' extends over all hadrons with pseudo-rapidity 
n < 3 in the Breit frame. The restriction to only the current hemisphere (77 < 0), as for the 2-jet shapes, would be 
too restrictive because of the extended phase space available for three partons. For reasons discussed in only 
events with p t ~ Q should be selected, which is accomplished by a cut on the (2+l)-jet resolution j/2 defined by the 
k t clustering algorithm: 0.1 < j/2 < 2.5. 

The azimuthal correlation between the hadrons labelled h and i is defined as 

X = 7T - \(f) h ~ (j>i\ , 
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where the observable is constructed by summing over all hadron pairs with a weight 

PthPu 



w 



Q 2 



The azimuth in the Breit frame of hadron h is denoted by 0/, . 

Predictions of mean values of 2-jet event shapes are available up to next-to-leading order (NLO) in the strong 
coupling, together with power corrections. Mean values of 2-jet rates for the JADE and the kt algorithm have 
been compared to NLO calculations, however the power corrections for this observables are not completely known. 
Soft gluon resummations to next-to-leading log precision (NLL) have been performed for distributions of 2-jet and 
3-jet event shapes. These have been matched to NLO fixed order predictions for the 2-jet event shapes only, power 
corrections are known for both. Distributions of jet rates can only be compared to NLO alone, as no matched soft 
gluon resummation nor power corrections are available yet. 



3. Mean Values 

In the Dokshitzer and Webber approach, the mean value of an event shape variable is modified through non- 
perturbative effects by an additive constant 

(F) = (F) pQCD + a F V, (8) 

where ap is of order one and can be calculated perturbatively. The power correction term V is assumed to be 
universal for all event shape variables. It is proportional to 1/Q and evaluated to be 



V = T M 15 

3n U 



ao&r) - a s (Q) - ^ fin £ + £ + 1 ) oJ(Q) 



2ir \ m (3 



(9) 



where 0q = 11 — 2n//3, K = 67/6 — ir 2 /2 — 5n//9, and n/ = 5 is the number of active flavours. The so-called 
Milan factor M. ~ 0.95 in the MS scheme ensures the universality at the two- loop level This ansatz results in 
only one single non-perturbative parameter ao = fij 1 f^ 1 a e g(k) dfc, being the first moment of the effective coupling 
integrated over the low scale region up to the matching scale pi . ao has to be determined from data and is expected 
to be around 0.5. 

For the perturbative part of the calculation the NLO programs DISENT 10] and DISASTER++ [HI have been 
used. The uncertainty due to missing higher orders of the perturbative series is by convention estimated by a variation 
of the renormalisation scale by a factor of two {fi r = Q/2 and fi r — 2Q). 



3.1 . 2-jet Event Shapes 

The HI Collaboration measured the mean values of five event shapes as function of Q , shown in Fig. A 
steady decrease of the means with rising Q is observed and a good description of the data is obtained by a fit of 
DISASTER++ supplemented by power corrections 0,^1 ■ F° r comparison the fixed order contribution alone is also 
given, which demonstrates the importance of the power corrections, especially for the non-photon-axis variables t c , 
Pq and C-Parameter. Earlier predictions of the jet broadening [1J| were not able to describe the data because of a 
more involved perturbative non-perturbative interplay of this observable. 

Similar results are obtained by the ZEUS Collaboration 0] , shown in Fig. [2 Note the logarithmic abscissa and 
the differing labels for the event shapes used compared to HI where 1 — Tt corresponds to t c , 1 — T 7 to r, B 1 to 
B and M 2 to po- Surprisingly the fitted power correction for the thrust with respect to the boson axis (1 — T 7 ) 
turns out to be slightly negative. The fitted values of a s (mz) and ao are shown in Fig. [3]and0] The 1-sigma 
contours denote experimental errors alone, which are only at the few percent level. These errors arc mainly caused 
by imperfect knowledge of the electromagnetic energy scale of the calorimeter, which enters the boost to the Breit 
frame of reference. Not shown are the sizeable theoretical errors (dominated by scale uncertainties due to missing 
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0( b 50 100 

Q/GeV 

Figure 1: Mean values of event shape variables corrected to the hadron level as a function of the scale Q. The data, presented 
with statistical errors (inner bars) and total errors (outer bars), are compared with the results of NLO QCD fits including 
power corrections (PC). The dashed curves show the NLO QCD contribution to the fits fl^| . 



higher orders) of about ten percent. In both analyses ao is found to be universal at the 10% level, and the spread of 
the fitted a s happens to be larger than expected from the experimental errors. There are some differences observed 
with the results from both collaborations, e.g. the values of ao fitted by HI cluster around a slightly higher value 
than those from ZEUS. However, the overall pattern looks quite similar, with r having a large error and a huge 
correlation between the fitted parameters. Also the observables which make not use of the virtual boson axis (t c ,pq 
and C-parameter) prefer rather high values of a s (mz) (compared to the world mean), in contrast to the broadening 
B, which exhibits the lowest values in both analyses. 

In contrast to e + e~ annihilation, where the perturbative coefficients are constant, in DIS they depend on x accord- 
ing to the parton density functions and the accessible x-range at different values of Q. Variables defined with respect 
to the boson axis, like r are expected to show a stronger dependence on x, which is demonstrated in Fig. |SJ While 
the perturbative contributions are expected to depend on x, the power corrections must not, in order of the concept 
to work. A former analysis by the ZEUS Collaboration Q observed an x dependent power correction parameter ao 
for r, t c and even po, however this observation was not confirmed by more recent investigations. 
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Figure 2: The mean values of event-shape variables as a function of Q . The solid lines are the results of the fit to the data of 
the predictions of the sum of NLO pQCD calculations from DISASTERS — h and the power corrections. The dashed lines are 
the DISASTER++ contribution to the fit alone |3. 

3.2. Jet Rates 

Jet rates show in general smaller hadronisation corrections than event shapes. Power corrections oc (A/Q) 2p have 
been proposed for the JADE and the k t jet algorithms, with suggested values of p = 1/2 and p = 1, respectively. 
Reliable calculations of the coefficients a/ have not been presented yet. In an HI analysis 17] an attempt was made 
to fit this coefficients together with a s and ao, the resulting theory curves are shown in Fig. |SJ While data are 
described quite well by the fit, the hadronisation corrections turn out to be negative and the fitted values of a s are 
unphysical. Since the effect of hadronisation is so small for this class of variables, even more precise experimental 
data could help to resolve this issue. 

4. Distributions 

Compared to mean values, the study of the spectra of event shape variables offers several advantages. Firstly, the 
shape of the distributions is governed by QCD and thus offers more information available for fits. Secondly, when 
fitting it is possible to restrict the range used to an interval where the theory prediction is reliable. 

4.1 . 2-jet Event Shapes 

Towards low values of the shape variables there are terms due to soft gluon radiation, which become dominant, 
but are not included in the NLO prediction used for the mean values. A significant improvement of the description is 
obtained with a resummation of these soft gluons, matched to the NLO part. Such a calculation is available in DIS 
for the 2-jet variables r, B, t c , po and C-Parameter. Fig. [7|shows exemplarily that a large part of the central region 
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Figure 3: Extracted parameter values for (a 3 ,ao) from fits to the mean values of the shape variables. The vertical line and 
the shaded area indicate the world average of a 3 (mz) 



of a spectrum is nicely reproduced by the theory, even at a rather low scale of Q = 15 GeV, where hadronisation 
effects are not small. 

The hadronisation correction in the Dokshitzer and Webber approach results in a shift for the differential distri- 
butions 

1 da{F ) _ 1 dgPQ CD (F-a F -p) ^ 
otot dF otot dF 

with the same coefficient of and power correction V as for the mean values. This shift cannot be valid over the 
whole spectrum, at low values of F it may be applied only for F clfV ~ fJ>i/Q Moreover, at large values of 
F higher order corrections are substantial and the perturbative calculation is not reliable. In consequence for each 
analysis one has to decide on the interval used for fits, depending on the observable and Q. For the upper bounds 
the values given in |t| are used, while for the lower bounds of the fit intervals different methods are applied. In case 
of the jet broadening the picture is more complicated, as the shift is supplemented by a squeeze [l9| . 

The calculations shown in the following were performed using DISASTER++ for the fixed order part and the 
DISRESUM package Q containing the resummed contributions and the power corrections. In a recent publication 
[lif the HI Collaboration has measured the differential distributions for the five event shapes quoted above, shown 
together with a QCD fit in Fig. IHlandED The data cover a wide range of (Q) = 15 — 116 GeV. Except for the 
highest Q bins, the precision of the measurements is not statistically limited. For each variable the shape of the 
spectra changes considerably with increasing Q, becoming narrower and evolving towards low values. The results of 
the combined fit for ao and a s (mz) are displayed in Fig. 1101 The quality of the fits, expressed in terms of \ 2 P er 
degree of freedom, is found to be between 0.5 and 1.4. It was checked that the results are stable against the omission 
of data points at the edges of the fit intervals. For all event shape variables consistent values for a s (mz) and «o are 
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Figure 4: Results of fits to the mean values of r, B, po, t~c and the C-parameter in the (a s , ao 
correspond to \ 2 = Xmm + 1> including statistical and experimental systematic uncertainties 
line) and its uncertainty (shaded band) are taken from [lg. 
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Figure 5: Mean values of r (left) and C (rig ht) versus Q in four different bins of x calculated with DISENT. The lines connect 
the means belonging to the same x bin |l7j|. 
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Figure 6: Mean values of y/j and yut as a function of Q The error bars represent statistical and systematical uncertainties. 
Upper part: The full line corresponds to a fit of the pQCD calculation without power corrections. Lower part: The full 
line corresponds to a power correction fit according to the Dokshitzer- Webber approach. The dashed line shows the pQCD 
contribution of DISENT in these fits [l^l- 

found, with a maximum difference of about two standard deviations between t and C . The fitted results for the five 
observables appear to be grouped for those with (r, B) and without (t c , po and C) referring to the boson axis in 
their definition. A strong negative correlation between a s (mz) and cvo is observed for all variables. The values of 
the strong coupling a s (mz) are in good agreement with the world average (ltjj. shown for comparison as the shaded 
band. The non-perturbative parameter ao — 0.5 is confirmed to be universal within 10%. A combined analysis of all 
event shape variables yields 

a s {m z ) = 0.1198 ±0.0013 (exp) tomlt (theo) , 

a = 0.476 ± 0.008 (exp) 1°;°^ (theo) , 

where the theoretical error is derived from the renormalisation scale uncertainty. The total errors are dominated by 
this renormalisation scale uncertainty, which suggests that missing higher order terms in the perturbative calculation 
are important. 
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Figure 7: Spectra of the Jet Broadening for (Q) = 15 GeV. The prediction by NLOJETH — h and soft gluon rcsummation (left) 
needs to be supplemented by power corrections in order to describe the data (right). 



A less consistent picture with respect to the universality of cxq is provided by preliminary results from the ZEUS 
Collaboration, shown in Fig. 1111 While the fitted values of a s (mz) are found to agree with the world average, a 
rather low value of cto is obtained for the C-Parameter, not covered by the experimental errors. However, only for 
r and B the value of \ 2 per degree of freedom is around unity, for the remaining variables it is around five. So the 
question arises whether the boundaries of the fit were stretched beyond the region where the theory is predictive. It 
is surprising that almost no (or even negative) hadronisation was found for the thrust (wrt. the boson axis) in case 
of the ZEUS mean value analysis (ao < 0.3 in Fig. [2J, while the differential distributions extracted from the same 
data to not show this property (with ao > 0.5). 

The event shape spectra used for the fit appear to be consistent with those from HI, but there are some differences 
in the treatment of the systematic errors. Moreover, there is some freedom when choosing options in the matching of 
the fixed order and the resummed part of the calculation and in deciding on proper fit intervals. Still, both analyses 
agree with respect to several findings: a s (mz) is compatible with the world mean (which is not the case for the event 
shape mean values), a s (mz) and cto have a negative correlation coefficient and the best description (quantified by 
% 2 /d.o.f.) is obtained for the boson axis variables. In order to draw more decisive conclusions about the validity 
of power corrections a consistent handling of the theory parameters, fit intervals and systematic errors would be 
desirable. Also, reduced theory uncertainties would help a lot, e.g. by a refined combination procedure [2flj | or even 
a NNLO calculation. 

In case of the HI analysis, the assumption of universal power corrections is assumed to be valid, and a fit as a 
function of Q, according to the renormalisation group equation is performed, with the resulting scale dependence of 
a s (Q) shown in Fig. 1121 



4.2. Jet Rates 

The longitudinally invariant k t algorithm is the de-facto standard for jet analyses in DIS. For this algorithm in its 
exclusive mode, the HI Collaboration determined distributions of 2-jet and 3-jet rates as a function of Q, shown 
in Fig. El Since the power corrections to this observables are not yet calculated, the data are compared to a NLO 
calculation, corrected for hadronisation effects by the Lund string model as implemented in the RAPGAP 2^1 event 
generator. At scales Q > 25 GeV a good description is provided for both jet rates, only at lower values deviations 
are found, which might be due to insufficient treatment of the hadronisation. However, as for the 2-jet event shapes, 
a soft gluon resummation matched to the fixed order part, might be important. Unfortunately, such a calculation is 
not available yet. 
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4.3. 3-jet Event Shapes 

Compared to 2-jet event shape, the 3-jet variables differ in that they are sensitive to large angle emissions between 
hard jets. Also it has been proposed to use such variables in a hadron-hadron collision environment. Preliminary 
measurements from the HI 21] and ZEUS Collaborations are shown in Fig. I14landll5l The theory predictions 
are not yet on the level of precision as for the 2-jet variables. A rather good description is obtained for the azimuthal 
correlation by the RAPGAP event generator, which relies on parton showers and the Lund string model. The out- 
of-plane momentum is compared to a LO+NLL calculation with power corrections, which does not a too bad job at 
high Q 2 . Some shift with respect to the data is seen, but one should keep in mind that the fixed order part of the 
calculation is only at leading order. As DIS data of good precision is available for both observables, the calculation 
at NLO+NLL precision would be highly desirable in order to draw further conclusions about power corrections. 

5. Summary 

Concerning the mean values of 2-jet event shapes overall consistent findings are reported by HI and ZEUS, with 
support for universal power corrections. A rather large spread in the fitted values of a s (mz) and the large scale 
uncertainties of the fitted results suggest that the currently available fixed order calculations at NLO are not sufficient 
for this application. 

For distributions of 2-jet event shapes soft gluon resummations are being used. In a recent HI publication the fits 
to the spectra resulted in compatible values for ao and the strong coupling among the five studied observables. In 
addition the obtained values of a s (mz) were found to be consistent with the world mean. A somewhat larger spread 
in ao is reported by a ZEUS preliminary result, but again the fitted a s [mz) are compatible with the world mean. 
Discrepancies in the findings of the two HERA experiments can be at least partly attributed to uncertainties due to 
determination of the fit intervals and due to differing matching prescriptions (between the fixed order the resummed 
calculations), as the measured data appear to be compatible with each other. Within the HI analysis it was shown 
that the results were not sensitive to changes in the fit interval. 

The topic of distributions from 3-jet event shapes is experimentally more difficult, because of reduced statistics and 
declined detector resolution. However, as good DIS data are on the market for two of these variables, a perturbative 
NLO+NLL could, when available, allow for a stringent test of power corrections beyond 2-jet variables. In case of 
the jet rates, besides a NLO+NLL calculation more insight into the power correction coefficients is eligible to draw 
further conclusions. 

On the other hand, the observed universality of ao for both distributions and mean values of the 2-jet event shape 
variables gives already strong support for the concept of power corrections. 

NB: After this workshop the ZEUS Collaboration made available new results on event shapes j^. The main 
characteristics of the fitted results have not changed considerably compared to the preliminary results shown in this 
article. 
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Figure 8: Normalised event shape distributions corrected to the hadron level for re, r and B. The data are presented with 
statistical errors (inner bars) and total errors (outer bars). The measurements are compared with fits based on a NLO QCD 
calculation including resummation (NLL) and supplemented by power corrections (PC). The fit results are shown as solid lines 
and are extended as dashed lines to those data points which are not included in the QCD fit fl^l . 
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Figure 9: Normalised event shape distributions corrected to the hadron level for po and the C-parameter. For details see the 
caption of Fig. |S] 
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Figure 10: Fit results to the differential distributions of r, B, po, re and the C-parameter in the (a a ,ao) plane lid. The 
la contours correspond to x 2 = Xmin + 1> including statistical and experimental systematic uncertainties. The value of a s 
(vertical line) and its uncertainty (shaded band) are taken from 16] . 
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Figure 11: Extracted parameter values for (a s ,ao) from fits to differential distributions of shape variables 1 1 el. The vertical 
line and the shaded area indicate the world average of n,(mz). 
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Figure 12: The strong coupling a s as a function of the scale Q. The individual fit results, shown as points with error bars, 
are obtained from fits to the differential distributions in tc, t, B, po and C within each Q bin. The errors represent the total 
experimental uncertainties. For each event shape observable a value of a 3 (mz) is indicated in the plot, determined from a fit 
to the a 3 (Q) results using the QCD renormalisation group equation. The corresponding fit curves are shown as full lines. The 
shaded bands represent the uncertainties on a s (Q) from renormalisation scale variations 
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Figure 13: The jet rate j/2 (left) and j/3 (right) |2l| compared to a calculation based on NLOJET++, where hadronisation 
corrections are determined with RAPGAP. 



R003 



15 



FRIF Workshop on First Principles Non-Perturbative QCD of Hadron Jets 



ZEUS 



ZEUS (prel.) 98-00 
LEPTO 

LO(0(a|))+NLL+PC 
100<Q 2 <500 GeV 2 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 14: K ou t distribution compared to LEPTO and the LO+NLL+PC calculation shown in bins of 
100 GeV 2 < Q 2 < 500 GeV 2 and 500 GeV 2 < Q 2 < 800 GeV 2 0. 
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Figure 15: Measured values of the azimuthal correlation |2l|. The data are compared with results from the RAPGAP event 
generator. 
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